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Abstract 

Third order WENO and CWENO reconstruction are widespread high 
order reconstruction techniques for numerical schemes for hyperbolic con¬ 
servation and balance laws. In their definition, there appears a small 
positive parameter, usually called e, that was originally introduced in or¬ 
der to avoid a division by zero on constant states, but whose value was 
later shown to affect the convergence properties of the schemes. Recently, 
two detailed studies of the role of this parameter, in the case of uniform 
meshes, were published. In this paper we extend their results to the case of 
finite volume schemes on non-uniform meshes, which is very important for 
h-adaptive schemes, showing the benefits of choosing e as a function of the 
local mesh size hj. In particular we show that choosing e = h~j or e = hj 
is beneficial for the error and convergence order, studying on several non- 
uniform grids the effect of this choice on the reconstruction error, on fully 
discrete schemes for the linear transport equation, on the stability of the 
numerical schemes. Finally we compare the different choices for t in the 
case of a well-balanced scheme for the Saint-Venant system for shallow 
water flows and in the case of an h-adaptive scheme for nonlinear sys¬ 
tems of conservation laws and show numerical tests for a two-dimensional 
generalisation of the CWENO reconstruction on locally adapted meshes. 

Keywords: WENO CWENO non-uniform mesh conservation and bal¬ 
ance laws high-order methods 
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1 Introduction 

Hyperbolic conservation and balance laws like 

u t + V • f(u) = g(u,x) (1) 
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describe a wealth of phenomena, from gas dynamics, to magnetohydrodynam¬ 
ics, to shallow water flows, that are of much practical interest in fields like the 
study of industrial processes, simulation-based prototyping, forecasting of nat¬ 
ural phenomena for the design of early-warning systems and passive defences 
against natural disasters, etc. Nonlinearities in the flux function /(it) give rise to 
very complex solutions, whose efficient numerical approximation often requires 
Adaptive Mesh Refinement (AMR) techniques, so that one can use fine meshes 
only in specific regions of the flow, e.g. close to boundaries, shocks, contact 
discontinuities. 

In the field of hyperbolic conservation laws, AMR has been achieved by 
superimposing patches of uniform grids, thus representing and time-advancing 
the solution at the same location on multiple meshes, taking care to maintain 
conservativity at mesh interfaces (3j. This approach undoubtedly has the ad¬ 
vantage that one can employ the very well-known and well-studied numerical 
schemes for discretization on uniform meshes. Alternatively, one may consider a 
single mesh, that can be locally refined, giving rise to an unstructured mesh. In 
this respect, it is important to notice that finite-volume methods can be easily 
formulated also on nonconforming meshes, like quad-tree ones, that can present 
hanging nodes. This approach, albeit facilitating the grid management and 
enforcement of conservativity, requires the development of new discretization 
techniques, suitable for irregular meshes. 

The fact that, in smooth regions of the flow, schemes of order three and above 
can compute accurate solutions on relatively coarse meshes is very important 
in AMR schemes, because a good adaptive strategy will be able to exploit this 
by leaving the mesh very coarse in smooth regions, and refining it only in a 
very small region around problematic features. For example, in [25] a second 
and a third order adaptive scheme with the same error indicator are compared, 
showing that the second order one has to refine in a much larger region, giving 
rise to much bigger meshes. 

Since a successful technique for mesh adaptivity in the case of hyperbolic 
systems is expected to continuously refine and coarsen the mesh at different 
locations in the computational domain, in our opinion it is important that dis¬ 
cretization techniques do not to rely on precomputed quantities related to the 
local mesh geometry (see e.g. [IT]). A good balance between the high accuracy 
and the small stencil requirement are third order accurate schemes. In this 
work we thus concentrate on two third order reconstructions that employ very 
small stencils, for which the needed metric information on the neighbours may 
be gathered efficiently at every timestep and extend to the non-uniform case a 
technique that can squeeze out from these small stencils a lot of extra accuracy 
with respect to the traditional form of the discretization. 

In order to fix the notation, we consider systems of balance laws with a 
geometric source term of the form 0 and we seek the solution on a domain f2, 
with suitable initial and boundary conditions. Although some two-dimensional 
applications will be shown in the numerical experiments at the end of the paper, 
the theory is developed in the one-dimensional case, in which the computational 
domain fl is an interval, discretized with cells f lj = ^+ 1 / 2 ], such that 
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U jQj = f 1 . The amplitude of each cell is hj = Xj + 1 / 2 — Xj_ i/ 2 , with the cell 
centre being Xj = (xj_ 1/2 + x j+ 1/2)/ 2 . 

We consider semidiscrete finite volume schemes and denote with 

u(t, x) dx 

the cell average of the exact solution in the cell Ij at time t and Uj(t) its 
numerical approximation. The semidiscrete numerical scheme can be written as 

~^Uj = (F j+ 1/2 - Fj_ 1/2) + Gj(U,x). ( 2 ) 

The numerical fluxes Fj + i/ 2 (t) should approximate f(u(t,Xj+ 1/2)) with suitable 
accuracy and are computed as a function of the so-called boundary extrapolated 
data , i.e. 

Fj+1/2 = F(U j+1/2 , U+ +1/a ), ( 3 ) 

where T is a consistent and monotone numerical flux, evaluated on two estimates 
of the solution at the cell interface U~^_ 1 j 2 , in turn obtained with a high order 
non oscillatory reconstruction. The function J- may be an upwind flux, (local) 
Lax-Friedriclis, an approximate or exact Riemann solver, etc. Finally, Gj is 
a consistently accurate, well balanced discretization of the cell average of the 
source term on the cell f lj, that in general requires the reconstruction of the 
point values of the solutions also at inner points in the cell (see e.g. PZOjh 

One of the most successful high order reconstructions is the Weighted Es¬ 
sentially Non-Oscillatory (WENO), whose first efficient implementation was de¬ 
scribed in [ 12 ]. It is based on the observation that one may recover the value of 
high order interpolating polynomials on a centred stencil as convex combination 
of the values of lower order ones that interpolate the function only in a subs¬ 
tencil (see [ 5 ] for a comprehensive development of this idea). In particular, let 
Pj,x(x) candidate polynomial in the j-th cell. Typically, Pj,x(x) will interpolate, 
in the sense of cell averages, the data U in a stencil A € A, where A is the set 
of stencil considered in the reconstruction. 

Next, depending on the particular reconstruction being performed, one has 
one or more sets of optimal weights C\ such that the linear combination C\Pj,\(x) 
has some desirable property, like coinciding with an higher order reconstruction 
at some point in the j-th cell. For example, the WENO reconstruction is char¬ 
acterised by two sets of optimal weights such that the combination of the 
evaluated polynomials AgA C±P\(xj ± |) matches the value of the central re¬ 
construction evaluated at the boundaries. The CWENO reconstruction instead 
has only one set of optimal weights and computes the linear combination of 
the candidate polynomials w hich can then be evaluated where 

needed. 

The weights of this convex combination are named linear weights but the 
actual reconstruction is computed by a modified set of weights, the so-called 
nonlinear weights , that are designed to be close to the linear ones in regions of 
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smoothness and close to zero if a discontinuity occurs in the corresponding sub¬ 
stencil. In order to achieve this, one considers suitable smoothness indicators, 
like those of [12J that defined 


DA = 2 

r> 1 


(h,r 


IQ 1 L 


P ^ X ) 


dx. 


(4) 


Note that these indicators are bounded even if the A-th stencil contain a dis¬ 
continuity, while they approach 0 if the solution is smooth in the stencil. In 
order to perform the reconstruction, one forms the nonlinear weights uj \, with 
the help of the regularity indicators as in 


WA 


C x 

(e + Ij,\) T 




Eje a s £ 


(5) 


and considers the linear combination The final result is a recon¬ 

struction that is close to some higher order one (dictated by the choice of the 
linear weights, often coinciding with the central one) in regions of smoothness 
and degrades smoothly to a one-sided non-oscillatory one close to discontinu¬ 
ities. Despite the absence of a formal proof of convergence for non- smooth 
solutions, this basic idea has been very successful and many details related to 
the very high order versions, the extension to higher dimensions, to non struc¬ 
tured meshes, the existence and positivity of the linear weights, together with 
dozens of applications, have been studied over the years. A comprehensive re¬ 
view with lots of useful references may be found in [25] . 

In order to obtain a fully discrete scheme, we apply a Strong Stability Pre¬ 
serving Runge-Kutta method (SSPRK) with Butcher’s tableau ( A , b ), obtaining 
the evolution equation for the cell averages 






(i) 

3 + 1/ 2 


- F 


W 


3 ~ 1/2 


+ At Y,b,G 


(i) 


( 6 ) 


Here Fj+ 1/2 = ^j+ 1 / 2 ) an d t ^ ie boundary extrapolated data U^’^ 2 

are computed from the stage values of the cell averages 
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k =1 


E 
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We point out that the spatial reconstruction procedures and the well-balanced 
quadratures for the source term must be applied for each stage value of the 
Runge-Kutta scheme. In this paper we consider a uniform timestep over the 
whole grid. A local timestep, keeping a fixed CFL number over the grid, can be 
enforced using techniques from mmm- The overall accuracy of the scheme 
is the smallest between the accuracy of the spatial reconstruction and that of 
scheme used for the time-evolution. 

With regards to the accuracy of the WENO reconstructions, (IQ) pointed out 
a previously unnoticed phenomenon of loss of accuracy close to local extrema 
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with non-vanishing third derivative and proposed a post processing (called map¬ 
ping) of the nonlinear weights aimed at closing the gap between linear and 
nonlinear weights for very smooth functions and recover the full accuracy also 
in this exceptional situation. Their original technique suitable for fifth order 
schemes, together with later improvements that extend it to the case of very 
high order schemes, is now known as WENO-M (see [7] and references therein). 
Another approach to solve the same problem, known as WENO-Z, consists in 
computing additional indicators, using the Jiang-Shu ones as building blocks, 
and it is effective from orders higher than 5 (see [61 and references therein). 

Recently, [Tl proposed a completely different technique. Upon observing 
that in the Jiang-Shu formulation m, the occurrence or not of the loss of 
accuracy is controlled by the relative size of the smoothness indicators and 
of the constant e, they proposed to set e = h 2 and showed the optimality 
of this choice in the case of uniform meshes, for WENO reconstructions of 
arbitrary order of accuracy. Indeed their technique allows to recover the correct 
asymptotic convergence rates even close to local extrema. Moreover, at the same 
time, it stabilises the convergence rates close to the asymptotic ones already 
for very coarse meshes, thus providing a reconstruction scheme that converges 
at the expected rate already on the meshes used in practice for computations 
and not only on very fine meshes. m extended this idea to the third order 
Compact WENO (CWEN03) scheme of [T7J, which was originally introduced 
in the context of central schemes since for the standard WEN03 scheme one 
cannot define linear weights that can yield third order accurate reconstructions 
at the cell centres. Of course the reconstruction is useful also in the case of 
non-central schemes for balance laws, when the reconstruction at the cell centre 
is needed in the well-balanced quadrature [UI1CI21- The CWENO reconstruction 
was extended to higher order and higher space dimensions mm, in the 
uniform grid setting and to non-uniform meshes of quad-tree type [231 . where 
the CWENO approach is very convenient since it is based on linear weights 
that do not depend on the relative size of the neighbours. In this last paper, 
the choice e = h was found useful in numerical experiments. WENO-M and 
WENO-Z are targeted to very high order reconstructions on uniform meshes, 
but in applications requiring local mesh refinement and unstructured meshes, 
due to the variability of the disposition and size of the neighbours and the fact 
that the mesh may change at every timestep, reconstructions with very small 
stencils and order up to four are to be favoured with respect to those requiring 
the collection of information over large stencils. This paper is devoted to the 
analysis of the non-uniform mesh case and compares the choices e = h 2 and e = h 
in the WEN03 and CWEN03 schemes in one space dimension. Several aspects 
are considered, from the formal accuracy, to stability, to the experimental orders 
of convergence observed on both coarse and fine meshes. 

In Section fj2] and )J3] we analyse the effects of different choices of e in the 
case of WEN03 and, respectively, CWEN03 schemes on non-uniform meshes. 
In f|4] we present numerical simulations corroborating our findings, including 
also tests on balance laws and an h-adaptive scheme for conservation laws. In 
i|5]we draw some conclusions and illustrate perspectives for future work. 
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2 WEN03 


The third order WENO reconstruction [21] considers the two first-neighbours 
of the cell in which the reconstruction is to be computed. We thus consider a 
non-uniform grid with hj = h , hj-\ = fih and hj+\ = 7 h. The two candidate 
polynomials are the linear polynomials matching the cell average of the central 
cell and of one of the neighbours, namely 


p jA x ) = u j + CT i,-0 - x i ). 
Pj,Fl{x) = Uj + CTj, + (x ~Xj), 


where <Tj t ± are given by 


- = 2 


Uj Uj — 2 

(1 + P)h ' 


a j,+ 


o U j +1 u j 

~ (1+7)* ’ 


(7a) 

(7b) 


The central reconstruction is the second order polynomial whose cell averages 
in f lj and match u ? and Uj±i ■ respectively, and is given by 

P J OPT (a;) = a + b(x — Xj) + c(x — Xj) 2 , (8) 

where 

_ — c 7 2 
a — Uj ^ 2 1 i 

, _ (5+/ 3 ) g j,+ + (|+7)g J .- 

l + /3 + 7 

3 (Tj ? -|- (7^ — 

2 /i(l + /3 + 7) 

One may easily verify that 

P? FT ( x j+ 1/2) = C~j, L Pj,L{ x j+ 1/2) + C'yfl- p i,fl( a; j+i/2) 
for 

r+ = 1 r+ = 1 + /3 

J,z ' 1 + (3 + 7 ’ J,ii 1 + /3 + 7 

and similarly that the optimal weights for the reconstruction at Xj- 1/2 are 


r~ = 1 + 13 r~ = 1 

x L 1 + /3 + 7 ’ i’ R 1 + /3 + 7 

Note that, in the case of a uniform mesh, /3 = 7 = 1 and one recovers the 

weights computed in [ 21 ] and pQ. 

In order to investigate the accuracy of the reconstruction, we apply it to the 
cell averages of a regular function u(x), whose Taylor expansions are given by 

Uj = Uj + fju" + ^ + ^(h 6 ), 

u j+ 1 = Uj + §(7 + 1 )u'j + |j(47 2 + 67 + 3 )u'j + |g (7 + 1 )( 27 2 + 27 + l)u'" + 0(h 4 ), 

1 = rtj — ^ (/? + 1)^ + ^j(4/3 2 + 6/3 + 3)uJ — (/3 -I-1)(2/3 2 -I- 2/3 -I-1 )u'” -t- 0(h 4 ). 
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The accuracy of the reconstruction of smooth data depends on the rate of 
convergence to zero of the discrepancy between the optimal and the nonlinear 
weights. 

Lemma 1. If u is smooth enough in the stencil Qj, then 


= 


(Uj ~ Uj- 1) 2 = 


(1 + /3) 2V J 

= u'j 2 h 2 - ^(2/3 + 1 )u' j u”h 3 + 

o 
4 


I i,R = 


(1+7 ) 5 

' 2 u2 


(Uj +1 - Uj) = 


= u'j h~ + -(27 + 1 )u'jUjh 6 + 


^[3(2 £ 2 + 2/3+ 1 )u' jU ”' + (4/3 2 + 4/3 + 1 )uf)h A + 0{h 5 ), 
— [3(2y 2 + 2y + 1 )u'jUj + ( 47 2 + 47 + 1 )u'j 2 }h 4 + 0(h 5 ). 


Remark 2. When the two regularity indicators Ij l and Ij r are very close 
to each other, procedure ([5]) gives nonlinear weights approximately equal to the 
linear ones and thus the reconstruction polynomial corresponds to the central sec¬ 
ond order polynomial. This is highly desirable for the reconstruction of smooth 
data. On the other hand, when a jump discontinuity is present, the imbalance 
between the regularity indicators will bias the reconstruction polynomial towards 
the smoothest stencil (e.g. uij,L — 1 and u>j,R — 0 if the jump is between Xj and 
Xj+i), avoiding oscillations but lowering the accuracy by one order. Unfortu¬ 
nately this may happen also close to local extrema of a smooth function, where 
both indicators would be close to 0 , but a slight imbalance between them may be 
amplified by the procedure © for the computation of nonlinear weights. 

Intuitively, ifu'(xj + 1 / 2 ) = Ij,R = (uj+i — Uj ) 2 would be much smaller 
than Ij } L = (Uj — Uj- 1) 2 and the nonlinear weights will differ quite a lot from 
the optimal ones, unless e is chosen such that it dominates both Ij^R and Ij t L 
in the denominators. Thus a value of e that does not depend on h will give 
rise to reconstructions that may change their order of convergence close to local 
extrema, depending on whether Ij,L/R e or not. The following computations 
will make this more precise and, for the above mentioned argument, we will 
concentrate on the case of non-constant h-dependent e. 


Theorem 3. Let u be a smooth function on a stencil {fb_i, Lli, with 

hj-\ = f3h, hj = h and hj +1 = 7 h. Let ojf L , uj 1 R be the nonlinear weights 
defined in (|5j) for the cell f lj. Then, we have 
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<L = C+ L 

E(Ev i+w “ >;v+o( H 

W 7,L = C i,L 

{ 1+ 3{l ^?Y‘ u ‘ hi+0(kt+1) \ 

5 

A 

So 

II 

£ 

[ 1 -3(«-E J ) 7 “ ; “ ;V + 0( '*‘ +I) . 
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So 

II 

E 

So 

1 - 777^727(1 + P»" hk + °^ fc+1 ) 
3{e + pu J ) 

= 0, k = 2 if e = 

eh , p = fc = 1 if e = eh 2 . 



Proof. Let 


M A E? - 


C i± 


j ’ R « ■ l,nr 

with r > 2. Then, following [I], we write 


C 


U i,L = 


0,L 


C 


j,L 


(e + I 0 ,l) t (e + Ij,rY 


1 + 


Ij,R - tj,L 


€ + I 


E 


i,L n 

s=0 


£ + Ij,R 
€ + I 


J,L 


By using the Taylor expansions ([9]), we obtain 

Yr-Yl = 2Y + H + 1) k + 0{hk+1) 
e + Y,L 3(e + pit' ) 

where k = 2 , p = 0 in the case e = eh and k = p = 1 if e = e/t 2 , thus 


E 

s —0 


£ + !>+ 

€ + Ij 


0,L 


l+s 2f p+ X u >> k+oihk+1) 

3(e + pit ) 


= E 

s =0 

r(r-l) 7 + /3 + 1 , „ k k+1 

= T + -o- ^—^ u i u 3 h +0{h + ). 

1 3(e + pit ■ ) 


Now, substituting in the expression (|12|, we have 

1 + 2 


= 




(e + -fj'.i ?) 7 


7 + P + 1 //, k 


^Tu'j u''h K + 0(h k+1 ) 


and 


w ,\£+ w jVR 


— ( £ + -fj.ii ) 7 


1 - 


3(e+pit' ) 

2 rip 


3 (e + pu'j ) 


J 2 -u' j U , 'h k + 0(h k+1 ) 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 


(14) 


with ip = 7 on the right of the cell, ip = 1 + ft on the left. Then, by using ( [IT] ), 
(13) and (141, with a simple computation, we obtain the desired results. □ 







































Remark 4. The result above establishes the Taylor expansion of ui\ —C\ up to 
the term that is enough for the proof of the third order accuracy of the recon¬ 
struction in Theorem [5| However, it is also interesting to observe the decay rate 
of uj\ — C\ close to a local extrema, i.e. to study the term 0(h k+1 ) in equation 
(10). Examining the proof above (see also m for the finite difference case), one 
can easily see that the accuracy order of uj\ — C\ is the same as that of j ( L . 
From <© we have that 

I J , R -I J ,L = l(T+P+l)u' J u''h 3 + ^-^( 1 +P + l)(3u' J u''' + 2uf)h 4 + 0(h 5 ) 

and Ij K = 0(h 2s+2 ), for K = L,R, where s = 0 if u)(xj) ^ 0 and s = 1 if 
u'jixjj = 0 £ 0 . 

So, for s = 0, we have 


0(h 3 ) 


Ij,R - Ij,L _ _ 

eh" + I jtL ~ eh" + (D{h 2 ) 


= 0{h 3 ~ v ), v = 0,1,2, 


with the case e = 10 ~ 6 corresponding to v = 0. On the other hand, for s = 1 
one finds 

Ij,R — Ij,L __ 0(h 4S| 

eh" + Ij,L 
in general, but 


eh" + 0(h 4 ) 

Ij,R ~ Ij,L . Q(h 5 ) 

eh" + I jL eh" + 0(h A ) 


= 0(h 4 -'), v = 0,1,2. 


= 0(h 5 -'), v = 0,1,2. 


in the case ft = 7 , which always occours on uniform grids. 
We denote with 


£-1/2 = Uj,L P jA X j-l/2) + u j,R P jA X 3-l/2)> 
l j+ 1/2 = u} t,L P iT( x j+ 1 / 2 ) + Wj' R Pj,R{Xj + 1 / 2 ), 


(15) 


the reconstructed values at the left and right boundary of cell flj and we prove 
a result on the accuracy of these boundary extrapolated values. 

Theorem 5. Let u be a smooth function in the stencil {f lj—i, Llj, STy+ 1 }. Then, 

u {. x j+ 1 / 2 ) = U j+l /2 0){h ) 
u ( x j-l/2) = u t-l/2 + 0( h3 )- 

Proof. We can compute the reconstruction error in the right side of the cell as 


l ( x j+ 1 / 2 ) - u j+i /2 ~ u ( x j+ 1 / 2 ) 


-p; 


OPT { 


( X j+1 /2) + P j ( x j+1 /2) - U j+l/ 2 - 


0(/t 3 ) 
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Indeed, the Taylor expansion gives 


Pr (x j+1/2 ) = u{x j+1/2 ) + ^(/3 + l)u'!'h* + 0(/r 4 ). 
Then, recalling that C+ R + C+ L = u>+ R + ui+ L = 1, and computing 

Pj,ni x j+ 1 / 2 ) = u(Xj + 1 / 2 ) + jr Ujh 2 + 0(h 3 ), 

fl + 1 //,2 


(16) 


(17) 


Pyi(a; j+1/ 2 ) = u(x j+1/2 ) - ——Ujh + 0(h 3 ), 


we estimate 


p 


OPT 


( Xj+1/2) U j+ 1/2 


(Cyi? _ - p 1,-r( ;c j+i/2) + (pf,L ~ -fj'P^i+i/s) - 

{pt,R _ w yi?) (-^'^(^j+i/a) ~ u ( x j+i/2)) + [pt, L ~ w ip) ( P h L ( x j+ 1/2) ~ ^(^+1/2))^ 


—v— 

0(h k ) 


0(h 2 


-V- 

0(h k ) 


- 

0(h 2 ) 


Using ( [To] ) in the above formula we have that k = 1 in the case of e = eh 2 
and k = 2 when e = eh. Thus, with both choices, the reconstruction error is of 
order 0(h 3 ). In a similar way, we can compute the reconstruction error in the 
left side of the cell f 1j. □ 

We note that, when e = eh 2 , all terms in the reconstruction error are of 
order 0(h 3 ), while, for e = eh, the interpolation error of P- PT dominates and 
the other terms are of higher order. This remark on the balance between the 
two sources of error is what makes p] state that the choice e = eh 2 is optimal 
and the computations above extend their result to non-uniform meshes. 

Next, we want to stress an important difference between finite difference 
and finite volume schemes. In fact, the accuracy of finite difference schemes 
depends on the accuracy of the derivative 1 - 1/2 — Fj- 1 / 2 ) which should 
be close to d x f(u{xj)), while semi-discrete finite volume schemes have their 
accuracy controlled by the discrepancy between P ?+1 /2 and f(u(xj + 1 / 2 )), see 
|161 §17]. In the following, we remark that, on uniform meshes, we have the 
correct accuracy of the derivative in the case of the linear transport equation, 
while in the numerical tests we will show that this may not maintained on non- 
uniform meshes, without affecting the convergence order of the fully discrete 
finite volume scheme. 

Remark 6. We have that in general 

(-1/2 H'j—1/2 1 / 2 ) u(Xj~ 1 / 2 ) + 0(h ), 

but it can be proven that there is an increased accuracy in the case of a uniform 
grid. In fact 


u j+ 1/2 u j- 1/2 


= (Pf^Xj+v 2) - PfJFixj-w)) 

+ (Pj- 1 (^-1/2) — U j_ 1/2) ~ (yPj ( X j+ 1/2) — U j + 1/2) 
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and for a uniform grid 

P ° PT ( x j+ 1 / 2 ) - P j 0 - P i T ( x j-i/ 2 ) = u(x j+1/2 ) - u(xj_ 1/2 ) + 0(h 4 ), (18) 

which raises the question if also the other terms can match this accuracy. 


( P j°- PT { x j-i/2) - V1/2) " { P ° PT ( x i+ 1/2) ~ ^+1/2) = 

— ~ W+fl) -fj>fl( a: j+l/ 2 ) “ (^ti “ -^^(^+1/2) 

and introducing the Taylor expansion of the approximation error of the linear 
polynomials the above expression becomes 

(c^_ 1R - vf_ UR ) (u( x j- 1 / 2 ) + \u”h 2 + 0(h 3 )) 

+ {c^_ 1L - (u(xj_ 1 / 2 ) - \ u 'jh 2 + 0{h 3 )) 

~ {pt,R ~ w yit) { u ( x j+ 1 / 2 ) + l u 'j h ~ + 0(h 3 )) 

- (CpL - Vj, l) («(*J+l/2) - + C’(^ 3 )) ■ 

Finally, collecting terms for each power of h and recalling that C^ R + C+ L = 
1 = u}j~ R + w+k, we find that the the approximation error is of order <D(h 4 ) if 

(“U ~ Hl) - (»U,R - Hit) = 0(h 2 ). ( 19 ) 


For e = eh the difference between linear and nonlinear weights are of order 
<D(h 2 ), see ( | 10 [ ) and thus the condition is trivially satisfied; furthermore, direct 
computation with e = eh 2 shows that condition (19) is satisfied also in this case. 

Note that condition (19) is analogous to the one found in m equation (29b)] 
for the fifth, order WENO reconstruction for finite difference schemes; there it 
was not found useful to design an optimal form of the nonlinear weights, but in 
this lower order case it is always satisfied. 


3 CWEN03 

The third order Compact WENO reconstruction (CWEN03) was introduced 
in na in the context of central schemes, originally in order to overcome the 
impossibility of finding linear weights for the reconstruction at cell centre using 
the WEN03 approach. It is of course useful also in non-staggered schemes, for 
example if a third order accurate reconstruction is needed at the centre of the 
cell, like in the well-balanced quadratures of [201 [22]. The main idea is to choose 
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( 20 ) 


linear weights Cj t L, Cj ,R G (0, l) such that Cj : L H - Cj,R ^ 1 and. to define 


Cj,0 = (1 - Cj )L - C j>R ), 

p j, oO) = (p° PT (x) - C jiL P jiL (x) - Cj t RP jtR {x))lCjfl, 

so that Pj- )PT (a;) = Cj > 0 Pj i 0 (x) + C^ L P^ L (x) + C jt RPj,n{x) at every point in 
the cell. Unlike in the WEN03 setting, here the values of Cj t L and Cj,R do not 
have to satisfy any accuracy-related requirement and thus can be in principle 
chosen arbitrarily; in particular we can employ the same linear weights, inde¬ 
pendently on the reconstruction point and on the local mesh geometry in the 
neighbourhood of the cell. 

With the choices Cj^ = Cjr = |, Cj t o = \ of [T7j one has that 
Pj,o{x) = Uj — ^ h 2 + (26 — tTj,+ + °' J ’~ )(x — Xj) + 2 c(x — Xj) 2 . 
Lemma 7. For u sufficiently smooth in the stencil flj, we have 


4,0 ~ 


+ 



a j,+ + g j, _|_ ^ c 2^4 _ u ,2 k 2 + |(/3 - 7 )u'ju'jh 3 + 
^j(0 - 7) 2 ^ n" 2 - ^ (/3 2 + 7 2 - 4/37 - /3 - 7 - 1 ) uju '" 


h 4 + 0(h 5 ). 


Theorem 8 . Let u he a smooth function on a stencil {Llj-i, Llj, with 

hj_ 1 = fih, hj = h and hj+\ = 7 h. We denote by uij^, Wj, R , uijfi the nonlinear 
weight defined in ([5]) in the cell f lj. Then , 


w jT - C j,L 


1 + _, 2/3 + 1 /2 Tu' jU ''h k + 0{h k+1 ) 


3 (e + pu'i ) 


Uj,R = C 


w j, 0 


j,R 


1 - 


27+1 


^T«'«"/i fc + 0 (/| t+1 ) 


fc+l'l 


: C 


j, 0 


1 + 


3(e + pu'. ) 
7-/3 


^rti'«;'/i‘+o(/i w ) 


3(e + pub ) 

with p = 0 , k = 2 for e = eh, and p = k = 1 for e = e 6 2 . 


( 21 ) 


Proof. Following the same procedure employed to get (101 in the WEN03 case, 
we obtain the Taylor expansions for the nonlinear weights. In particular, for 
Wj t 0 , we find, 


Uj,0 - 


Cj, 0 

(e + 4 o) T 


Cj, 0 

(e + 4,4 T 



By using the Taylor expansions ([ 9 ]), we obtain 


Ij,R Ij, 0 

e + 4,0 


37 — /3+1 
3(e + pub 2 ) 


u' jU ”h k + 0(h k+1 ) 


( 22 ) 
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where k = 2, p = 0 if e = eh and k = p = 1 if e = eh 2 . Thus, 


E 

s—0 


e + fj'R 

e + Ij, 0 


= £ 1 + » ?'. 9 tv “>?'■* + °('» w ) 

s =o L 3 ( e +P u j) 

- r+ T(T D3-, ^ + 1 (AW) 

2 3(e + pu'- J 


Now, substituting in (22) 

G*,o 


w i,o - 


(e + 


1+ V-4.' S ^ T "i"?'‘ t+0(l ‘ W) 

3(e + pu', ) 


and 


1 


Wj,L + ^j,R + w i,0 


— (e + Ij,R) T 


1 - 


2 7 + l 


3 {e+pu'i ) 


7 2 -TU , j v!!h k + 0{h k+1 ) 


Then, recalling ( flTj ) , ( [l3| , with the different values for Cjj J and C^r, we obtain 
the results by substituting in 


w i,o - 


Uj, o 


“i,L 


Uj,R 


<^j,L + u j,R + w i,0 + Wj,.R + UJ^o ’ w j,L + + WjjO 


□ 


Remark 9. Using the same notation of Remark^ in this case we have that 


Ij,R - Ij ,o = g(3 7 - P + 1 )ujUj h 3 + 


1 

12 


, a ^, 0 Q , ^ / /// , 3 7 2 -/? 2 + 4 7 + 2/3 7 - 155 //2 ' 

( 7 -p)(3 7 -/3 + l)M j w j H- - - Uj 


0(h 5 ) 


and thus Ij r — Ij o = 0(fh 3+s ) independently of the grid employed. With the 
same argument of Remark ul using Ij o = 0(h 2s+2 ), one shows that, when 

u' j (x j )=0, C x -to x = 0(h^)- 

In the CWEN03 reconstruction, only one set of nonlinear weights is com¬ 
puted and then uj +1 , 2 and uf_ 1 , 2 are simply obtained by evaluation of Pj(x) = 
<jJj,LPj,L{x) + uij^RPj } B.{x) + ojjfiPjp(x) at the boundaries of the cell. 

Theorem 10. Let u be a smooth function in the stencil {ST,-_i, flj, flj+i}. Then, 

u(x j+1/2 ) = uJ +1/2 + 0{h 3 ), 
u{xj_ i/ 2 ) = u+_ 1/2 + 0(h 3 ). 
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Proof. Similarly to the WEN03 case, from (161, we can compute the recon¬ 
struction error as 


u ( x j+ 1 / 2 ) - u j+ 1/2 — u(xj+ 1 / 2 ) - Pj >PT (x j+1 / 2 ) +P J OPT (a: : , + i/ 2 ) - Wj- +1 / 2 - 

0 (h 3 ) 


Then, recalling that Cy# + Cj t L + Cj ,0 = + w j,o = 1 ; from p~7| , and 

using the Taylor expansion of the quadratic polynomials 

Pj,o( x j+i/2) = u{xj + 1/2) H-^+ ^(^ 3 )) (23) 


we estimate 


Pj ( x j+ 1/2) M j+1/2 

(Cj'.-R _ -fj,-R( :E i+i/2)+(Cj',-L _ w i,i) -Pi,i( a: i+i/2)+(C'i,o - Vj,o) -Pyo (£3+1/2) = 

(Cyfl - Wy R ) (Pj, R (x j+1/2 ) - u(x j+ 1/2)) + (Cyz, - wyz,) (Pyi(^ + i/2) - m(^ + i /2 )) + 

v -V--V-' v -V-''-V-" 

0 (h k ) 0 (h 2 ) 0 (h k ) 0 (h 2 ) 

+ (Cj,0 — w yo) (-Pyo(£j+i/2) — ^(^+1/2)) ■ 

'- v --- v _ ' 

0 (h k ) 0 (h 2 ) 


Thanks to (21), in the formula above k = 1 in the case of e = eh 2 , while k = 2 
when e = eh. □ 


We note that, with both choices, the reconstruction error is of order 0(h 3 ), 
but, when e = eh 2 , all terms in the reconstruction error are of order 0(h 3 ), while 
for e = eh, the interpolation error of P J OI>T dominates and the other terms are of 
lower order. These computations extend the results of M to the non-uniform 
case. 

We conclude this section with the analogous of Remark [6] in the case of 
the CWEN03 reconstruction. As for the previous case, this result should be 
contrasted with the numerical evidence of §4.2| 

Remark 11. Proceeding as in Remark [fi| we obtain that, using a uniform grid, 

u J+i /2 ~ V- 1/2 = u ( x j+ 1 / 2 ) - u(xj_ 1/2 ) + 0(h 4 ). 

In the CWENOS approximation we have 


(P°- P 1 T ( X 3-1 /2) - V 1 / 2 ) - ( P ° FT (*3 +1 /2) - ^ + l/ 2 ) = 

(Cj- 1 ,R ~ -l,n) Pj-l,R{ x j-l/ 2 ) + {Cj-t,L ~ Pj-l,L{Xj-l/ 2 ) 

T ^j— 1,0) Pj—i,o(,Xj— 1/2) (P'j, 0 Mj, 0) Pj,oi.Xj+ 1/2) 

- (Cj,R - UJj, R ) Pj, R (x j+ 1/2) - (Cj'L - Uj,L) Pj,L{x j+ 1/2). 
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Using the Taylor expansion of the quadratic polynomials 


Pj-i,o(xj-i/2) — u(xj_ 1/2) + + 0 [h 3 ), 

Pj, o{ x j+X/l) = u{Xj + 1/2) + J 2 u 'jh 2 + @(h 3 ), 


collecting terms for each power of h, recalling (18) and that Cj t u + C 3i l + Cj ,0 = 
u)j t R+ojj t L+ojj t o = 1 , we find that the the approximation error is of order 0(h 4 ) 

if 


(&j,o — Vj- 1,0) + 2 [ — 2 u)j t L) — — 2 a; J -_i i i) ] = 0 {h 2 ). 


(24) 


-For e = eh the difference between linear and nonlinear weights are of order 
0(h 2 ), see ( | 21 [ > with k = 2 , so the condition is trivially satisfied; furthermore, 
in the case e = eh 2 , the condition (24) is satisfied by using (21) with k = 1. Note 
that equation (24) is analogous to (19) in the case of CWEN03 reconstruction. 


4 Numerical tests 


First we describe the meshes used in the numerical tests, on the reference interval 
[0,1]. Uniform grids with N cells, have cells of size h = 1/N and cell centres 
Xj + i / 2 = (j + \ )h. A quasi regular grid is obtained transforming the cell centres 
of a uniform grid with the map 


<p{x) = x + 0.1 * sin(107ra;)/5. 


The grid spacing in quasi-regular grids is asymptotically described by In 

this way, we obtain a grid with cell sizes hj such that 

(1 ~ f )w — — (1 + f )~k- 


Note that the ratio of consecutive cells approaches 1: 


h k+1 1 y"(£) 

h k N V'(w)' 


(25) 


In the notation of Sections [2] and [3] this means that, for increasing N, the 
parameters (3 and 7 approach the value 1 that characterise uniform meshes. 
We will observe that the numerical schemes on quasi-regular grids behave very 
much like on uniform ones. 

Next, we consider random grids that are obtained moving randomly the 
interfaces of a uniform grid. We start from a uniform grid with cells of size 
h = 1/N and then we consider grids with interfaces at 


Xj + l/2 — Xj+ 1/2 + 


where are random numbers uniformly distributed in [—0.5,0.5]. We have 


il < h j < 


5 J_ 
4 N * 
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We use this grid for the purpose of illustration even if, of course, one would not 
use such irregular grid in an application. 

Finally, observe that in binary-tree mesh refinement in d spatial dimensions, 
one starts with a uniform mesh and then, guided by some error indicator, each 
cell may be split (recursively) in 2 d equal parts, see [211 . The ratio of 

the sizes of adjacent cells is thus a positive power of 2 and does not approach 
1. In order to test the schemes on meshes that may be employed by such 
AMR techniques, we consider the a/3^5 meshes, that are composed by adjoining 
building blocks of 5 cells of size ah, j3h, h , yft,, Sh respectively. 

4.1 Comparing e’s 

The first set of tests is on the reconstruction error of WEN03 and CWEN03 on 
non-uniform meshes. An a/Sytj-grid is set up with cells of size h, 2 h, h, h/ 2, h/ 2 
with x = 0 at the centre of the middle cell. In all our tests, we consider nonlinear 
weights ([5]) with r = 2. For decreasing values of h, cell averages of a function u(x) 
were set using the fifth order gaussian quadrature rule and the reconstruction 
at the right boundary of the middle cell was computed and compared with the 
exact values u(h/2). 

In Table [I] we show the reconstruction errors and convergence rates for 
WEN03 on non-uniform meshes, with different choices of e. It is clear that 
the non-uniformity of the mesh has no influence on the convergence rates. In 
fact, as in the uniform grid tests of [I], an irregular convergence rate appears for 
constant e when u' vanishes in the central cell, but a regular convergence history 
can be recovered by employing an /i-dependent e. We remark that e = h gives 
slightly lower errors than e = h 2 . We also point out that repeating the test of 
Table [l] for a function such that u'(h/2) = 0 , gives analogous results, indicating 
that convergence may be degraded whenever v! vanishes in the reconstruction 
stencil. 

In Table [2] we show the same tests for CWEN03; we obtained analogous 
results and the remarks about the comparison with the uniform grid case could 
be repeated. 

In Figure [l] we show the distance |wa — C\ | observed when reconstructing 
u(x) = x 3 + cos x on a non-uniform grid of type a, /3, 7 , 5. As expected, in both 
the WEN03 and CWEN03 cases, the choice e = 10“ 3 ° does not give weights 
converging to their optimal values, the choice e = 10 -6 behaves similarly on 
coarse meshes and changes to a convergent regime at about h = 2 x 10 -3 . On 
the other hand, the choices e = h and e = h 2 give rise to a more regular con¬ 
vergence histories, with the former giving lower discrepancies between nonlinear 
and optimal weights. 

Furthermore, Tables [3] and [4] report the asymptotic behaviour of the nonlin¬ 
ear weights for the WEN03 and CWEN03 recontruction, comparing the case 
u'j 7 ^ 0 and the case with a smooth extremum in the reconstruction stencil. The 
behaviours follow the results of Remarks 0] and [9] 

Finally, in order to investigate more deeply the behaviour of the reconstruc¬ 
tion procedure, we denote by , Uj, %+i) the map from the three cell av- 
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Reconstruction error for u(x) = e x 


h 

e = 10- 

cot 

e = 10- 

u 

e = h 

e = h 

2 


error 

rate 

error 

rate 

error 

rate 

error 

rate 

5.00e-02 

1.16e-05 


1.16e-05 


2.31e-06 


4.70e-06 


2.50e-02 

1.43e-06 

3.02 

1.42e-06 

3.03 

3.08e-07 

2.91 

5.65e-07 

3.06 

1.25e-02 

1.78e-07 

3.01 

1.72e-07 

3.04 

3.96e-08 

2.96 

6.92e-08 

3.03 

6.25e-03 

2.21e-08 

3.00 

1.96e-08 

3.13 

5.02e-09 

2.98 

8.56e-09 

3.01 

3.12e-03 

2.76e-09 

3.00 

1.78e-09 

3.47 

6.32e-10 

2.99 

1.07e-09 

3.01 

1.56e-03 

3.45e-10 

3.00 

8.15e-ll 

4.45 

7.92e-ll 

3.00 

1.33e-10 

3.00 

7.81e-04 

4.31e-ll 

3.00 

2.92e-12 

4.80 

9.92e-12 

3.00 

1 .66e-ll 

3.00 

3.90e-04 

5.38e-12 

3.00 

9.99e-13 

1.55 

1.24e-12 

3.00 

2.07e-12 

3.00 

1.95e-04 

6.73e-13 

3.00 

1.48e-13 

2.78 

1.55e-13 

3.00 

2.59e-13 

3.00 

9.76e-05 

8.39e-14 

3.00 

1.91e-14 

2.95 

1.95e-14 

3.00 

3.22e-14 

3.00 


Reconstruction error for 

u{x) 

= cos (27ra:) + x 3 



h 

e = 10- 

TsU 

e = 10- 

6 

e = h 

e = h 

2 


error 

rate 

error 

rate 

error 

rate 

error 

rate 

5.00e-02 

7.91e-03 


7.91e-03 


7.61e-04 


6.79e-03 


2.50e-02 

2.00e-03 

1.99 

1.99e-03 

1.99 

3.12e-05 

4.61 

1.06e-03 

2.68 

1.25e-02 

5.01e-04 

2.00 

4.75e-04 

2.07 

1.41e-06 

4.47 

9.72e-05 

3.45 

6.25e-03 

1.25e-04 

2.00 

4.91e-05 

3.28 

8.19e-08 

4.10 

6.77e-06 

3.84 

3.12e-03 

3.13e-05 

2.00 

1.04e-06 

5.55 

6.35e-09 

3.69 

4.36e-07 

3.96 

1.56e-03 

7.84e-06 

2.00 

1.71e-08 

5.93 

6.14e-10 

3.37 

2.77e-08 

3.98 

7.81e-04 

1.96e-06 

2.00 

3.26e-10 

5.72 

6.75e-ll 

3.19 

1.76e-09 

3.97 

3.90e-04 

4.90e-07 

2.00 

1 .20e-ll 

4.77 

7.92e-12 

3.09 

1.14e-10 

3.95 

1.95e-04 

1.22e-07 

2.00 

1 .02e-12 

3.55 

9.60e-13 

3.04 

7.59e-12 

3.91 

9.76e-05 

3.06e-08 

2.00 

1.19e-13 

3.10 

1.18e-13 

3.01 

5.33e-13 

3.83 


Table 1: WEN03. Reconstruction errors at x = 0 + h/2 for a grid of five cells 
of size h, 2h, h, h/2, h/2 with x = 0 in the centre of the middle cell. In the first 
test v! ^ 0 in the reconstruction stencil, while u , (0) = 0 in the second case. 
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Reconstruction error for u(x) = e x 


h 

e = 10- 

error 

cot 

rate 

e = 10- 

error 

rate 

e = i 

error 

h 

rate 

e = h 

error 

/2 

rate 

5.00e-02 

4.60e-06 


4.59e-06 


2.50e-06 


1.05e-06 


2.50e-02 

5.58e-07 

3.04 

5.53e-07 

3.05 

3.19e-07 

2.97 

1.19e-07 

3.13 

1.25e-02 

6.88e-08 

3.02 

6.61e-08 

3.06 

4.03e-08 

2.99 

1.42e-08 

3.07 

6.25e-03 

8.54e-09 

3.01 

7.28e-09 

3.18 

5.06e-09 

2.99 

1.74e-09 

3.03 

3.12e-03 

1.06e-09 

3.01 

5.70e-10 

3.67 

6.34e-10 

3.00 

2.15e-10 

3.02 

1.56e-03 

1.33e-10 

3.00 

9.70e-13 

9.20 

7.94e-ll 

3.00 

2.67e-ll 

3.01 

7.81e-04 

1.66e-ll 

3.00 

6.43e-12 

-2.73 

9.93e-12 

3.00 

3.32e-12 

3.00 

3.90e-04 

2.07e-12 

3.00 

1.12e-12 

2.52 

1.24e-12 

3.00 

4.15e-13 

3.00 

1.95e-04 

2.59e-13 

3.00 

1.52e-13 

2.88 

1.55e-13 

3.00 

5.15e-14 

3.00 

9.76e-05 

3.22e-14 

3.00 

1.95e-14 

2.96 

1.95e-14 

3.00 

6.44e-15 

3.00 


Reconstruction error for 

u{x) - 

= cos (27ra;) + x 3 



h 

e = 10- 

TsU 

e = 10- 

-6 

e = h 


e = h : 

> 


error 

rate 

error 

rate 

error 

rate 

error 

rate 

5.00e-02 

7.85e-03 


7.85e-03 


4.81e-04 


6.38e-03 


2.50e-02 

1.98e-03 

1.99 

1.98e-03 

1.99 

2.05e-05 

4.56 

8.49e-04 

2.91 

1.25e-02 

4.97e-04 

2.00 

4.64e-04 

2.09 

1.07e-06 

4.27 

6.06e-05 

3.81 

6.25e-03 

1.24e-04 

2.00 

3.58e-05 

3.70 

7.11e-08 

3.91 

3.65e-06 

4.05 

3.12e-03 

3.11e-05 

2.00 

5.48e-07 

6.03 

6.01e-09 

3.56 

2.25e-07 

4.02 

1.56e-03 

7.78e-06 

2.00 

8.89e-09 

5.94 

6.04e-10 

3.31 

1.42e-08 

3.98 

7.81e-04 

1.94e-06 

2.00 

1.96e-10 

5.50 

6.72e-ll 

3.17 

9.16e-10 

3.95 

3.90e-04 

4.86e-07 

2.00 

9.93e-12 

4.31 

7.92e-12 

3.09 

6.10e-ll 

3.91 

1.95e-04 

1.22e-07 

2.00 

9.91e-13 

3.32 

9.60e-13 

3.04 

4.28e-12 

3.83 

9.76e-05 

3.04e-08 

2.00 

1.19e-13 

3.06 

1.18e-13 

3.01 

3.25e-13 

3.72 


Table 2: CWEN03. Reconstruction errors at x = 0 + h/2 for a grid of five cells 
of size h, 2h, h, h/2, h/2 with x = 0 in the centre of the middle cell. In the first 
test u' ^ 0 in the reconstruction stencil, while u , (0) = 0 in the second case. 
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u 

(x) = x 3 + cos(27ra;) 



e = 10- 

a 

e = h 


e = h 2 


h 

Cr ~ 

rate 

Cr - <^R 

rate 

Cr - <^R 

rate 

5.00e-02 

1.38e-01 


2.27e-02 


1.67e-01 


2.50e-02 

1.38e-01 

0.00 

3.37e-03 

2.75 

7.45e-02 

0.65 

1.25e-02 

1.37e-01 

0.02 

4.33e-04 

2.96 

2.89e-02 

1.37 

6.25e-03 

1.06e-01 

0.36 

5.44e-05 

2.99 

8.29e-03 

1.80 

3.12e-03 

1.89e-02 

2.49 

6.80e-06 

3.00 

2.15e-03 

1.95 

1.56e-03 

1.32e-03 

3.84 

8.50e-07 

3.00 

5.42e-04 

1.99 

7.81e-04 

8.29e-05 

3.99 

1.06e-07 

3.00 

1.36e-04 

2.00 

3.90e-04 

5.18e-06 

4.00 

1.33e-08 

3.00 

3.40e-05 

2.00 

1.95e-04 

3.24e-07 

4.00 

1.66e-09 

3.00 

8.49e-06 

2.00 

9.76e-05 

2.02e-08 

4.00 

2.07e-10 

3.00 

2.12e-06 

2.00 




u(x ) = e 

,X 




e = 10- 

-6 

e = h 


e = h z 


h 

Cr — UR 

rate 

Cr — <^R 

rate 

Cr — <^R 

rate 

5.00e-02 

3.08e-02 


1.33e-03 


1.47e-02 


2.50e-02 

1.48e-02 

1.06 

3.44e-04 

1.95 

7.22e-03 

1.02 

1.25e-02 

7.23e-03 

1.03 

8.75e-05 

1.97 

1.79e-03 

1.01 

6.25e-03 

3.51e-03 

1.04 

2.21e-05 

1.99 

8.94e-04 

1.00 

3.12e-03 

1.63e-03 

1.11 

5.55e-06 

1.99 

4.47e-04 

1.00 

1.56e-03 

6.34e-04 

1.36 

1.39e-06 

2.00 

2.23e-04 

1.00 

7.81e-04 

1.69e-04 

1.91 

3.48e-07 

2.00 

1.12e-04 

1.00 

3.90e-04 

2.95e-05 

2.52 

8.71e-08 

2.00 

5.58e-05 

1.00 

1.95e-04 

4.10e-06 

2.85 

2.18e-08 

2.00 

2.79e-05 

1.00 

9.76e-05 

5.27e-07 

2.99 

5.45e-09 

2.00 

1.39e-05 

1.00 


Table 3: WEN03. Distance and order of convergence of the optimal weights 
Cft and the nonlinear weights , as a function of the mesh width h, in the 
cell centered in x = 0. 
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u(x) = X s + cos( 27 rx) 

h 

e = 10- 


e = 

h 

e = h 1 



Co — wo 

rate 

Co — wo 

rate 

Co — w>o 

rate 

5.00e-02 

4.99e-01 


2.43e-01 


4.89e-01 


2.50e-02 

4.99e-01 

0.00 

4.55e-02 

2.42 

4.42e-01 

0.15 

1.25e-02 

4.98e-01 

0.00 

6.02e-03 

2.92 

2.82e-01 

0.65 

6.25e-03 

4.81e-01 

0.05 

7.59e-04 

2.99 

1.05e-01 

1.43 

3.12e-03 

2 .10e-01 

1.20 

9.50e-05 

3.00 

2.93e-02 

1.84 

1.56e-03 

1.81e-02 

3.53 

1.19e-05 

3.00 

7.53e-03 

1.96 

7.81e-04 

1.16e-03 

3.97 

1.48e-06 

3.00 

1.90e-03 

1.99 

3.90e-04 

7.24e-05 

4.00 

1.85e-07 

3.00 

4.75e-04 

2.00 

1.95e-04 

4.53e-06 

4.00 

2.32e-08 

3.00 

1.19e-04 

2.00 

9.76e-05 

2.83e-07 

4.00 

2.90e-09 

3.00 

2.97e-05 

2.00 




u(x ) = 

= e x 



h 

e = 10- 

-6 

e = 

h 

e = h z 



Co — wo 

rate 

Co — oj o 

rate 

Co — wo 

rate 

5.00e-02 

3.16e-02 


1.40e-03 


1.52e-02 


2.50e-02 

1.42e-02 

1.16 

3.32e-04 

2.08 

6.96e-03 

1.13 

1.25e-02 

6.64e-03 

1.10 

8.06e-05 

2.04 

3.30e-03 

1.07 

6.25e-03 

3.15e-03 

1.07 

1.98e-05 

2.02 

1.61e-03 

1.04 

3.12e-03 

1.44e-03 

1.13 

4.92e-06 

2.01 

7.93e-04 

1.02 

1.56e-03 

5.59e-04 

1.37 

1.22e-06 

2.01 

3.93e-04 

1.01 

7.81e-04 

1.49e-04 

1.91 

3.06e-07 

2.00 

1.96e-04 

1.01 

3.90e-04 

2.59e-05 

2.52 

7.64e-08 

2.00 

9.78e-05 

1.00 

1.95e-04 

4.59e-06 

2.85 

1.90e-08 

2.00 

4.89e-05 

1.00 

9.76e-05 

4.61e-07 

2.96 

4.77e-09 

2.00 

2.44e-05 

1.00 


Table 4: CWEN03. Distance and order of convergence of the optimal weights 
Cq and the nonlinear weights ujq , as a function of the mesh width h, in the cell 
centered in x = 0. 
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Figure 1: Distance between the optimal weights C\ and the nonlinear weights 
uj\ as a function of the mesh width h, for different choices of e. WEN03 (left) 
and CWEN03 (right) 


erages to the reconstructed boundary value Uj + i/ 2 and by V 2 {uj-\,Uj,Uj+\) = 
Pj )FT (Xj + 1/ 2 ) the reconstruction operator that employs the central parabola. 
The previous test confirmed that 1Z(uj-i,Uj, Uj+i) —>• V 2 (uj-i,Uj,Uj+i) for 
h — > 0. In Figure [2] we test the convergence \71Z —» V"p2- The plots show that 
for /i-dependent e, dlZ/duk converge quickly to &P 2 / duk for k = j — 1 ,j,j + 1. 
At the opposite hand, for e = 10 -30 we do not observe such a convergence (note 
that dlZ/duk -f* cfP 2 /<9ufc, k = j — 1, j,j + 1) and e = 10 -6 shows an hybrid 
behaviour that changes regime when h falls below a threshold. (These tests 
were performed on uniform meshes for u(x ) = x 3 + cosx). 

4.2 Numerical derivative and linear transport 

In this set of tests we investigate the effects of the choice of e in a numerical 
scheme for the linear transport equation. Recall that, for ut + u x = 0, when 
using upwind numerical fluxes, the semidiscrete scheme © boils down to 

$-11; =-^(U~ /0 - U7 , . 

df J hj V J+ 1 / 2 J-V2; 

This system of ODE is discretized with the third order, three stages SSPRK, 
see [9]. For these tests, a non uniform grid obtained repeating groups of four 
cells of size h, h/ 2, h/ 4, h/ 4 was generated. We point out that cells of size h 
have a neighbourhood of type a, (3,^,6 = 1/4, 1/4, 1/2, 1/4, cells of size h/2 
have a neighbourhood with a,/?, 7, S = 1/2, 2, 1/2, 1/2, while cells of size h/ 4 
have a, / 3 ,7, S = 4, 2, 1, 4 or a, f), 7, S = 2, 1, 4, 2. 

We recall that, if the numerical flux function T appearing in © is consistent 
and Lipschitz continuous, the reconstructions are third order accurate and the 
ODEs ([2]) are discretized using a third order accurate method, then the numer¬ 
ical finite volume scheme is third order accurate in time as well in space, see 

hei § 17 ]. 
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Figure 2: Difference between the components of 'S/TZ(uj-i,Uj,Uj+i) and of 
VP2 (uj-i, Uj , u j+i ) as a function of the uniform mesh size h for different choices 
of e. WEN03 (left) and CWEN03 (right). 


These tests compute also the spatial discretization error 


(u{x j+1/2 ) - m (^_ 1 / 2 ))// i j - (U j+1/2 - U j l/2 )/hj, 

which is the finite volume error analogue of the finite difference truncation error 
“'(*>■) - ("w - u s -1/2) /hj Studied in 1 ). 

We integrate Ut + u x = 0 on the domain [0,1] until t = 1, with peri¬ 
odic boundary conditions and the smooth initial datum uq(x) = sin(27ra — 
sin(27ra;)/27r). The maximum norm error for the numerical derivative of Uq(x) 
on this grid and the 1-norm error at final time in the linear transport test were 
recorded. The results for e = h 2 are shown in Table [5] (WEN03) and Table 
[6] (CWEN03). Uniform and quasi-regular grids showed the expected conver¬ 
gence rates ((251 and Remarks [6] and 11) in both the spatial discretization error 
and the linear transport test. Random and a/SyiJ-grids show irregular and de¬ 
graded convergence rates for the spatial discretization error, but the order of 
convergence is 3 in the linear transport test. 

Having shown that random and a/3yi5-grids are the most troublesome, but 
that in the linear transport test the theoretical order of convergence is easily 
reached with an /i-dependent e, next we compare the different combinations of 
WEN03 and CWEN03 with e = h 2 and e = h on linear transport test for 
smooth and discontinuous data. 

The error at final time in both the 1-norm and the maximum norm are 
reported in Figure [3] On all grid types (only random ones are shown) and for 
both norms, the choice e = 1C) -6 gives the biggest errors, e = h 2 is slightly 
better, while e = h yields errors that are lower by about a factor of 2. In all 
cases, the CWEN03 reconstruction yields slightly lower errors than the WEN03 
reconstruction. 

Next we repeat the previous test with the square wave initial datum uq(x) = 
X[i/ 2 ,i](x), where x denotes the characteristic function. At final time we com- 
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Maximum norm error on the numerical derivative on uq(x) 


N 

uniform 

error rate 

quasi¬ 

error 

unif 

rate 

random 

error rate 

[1,1/2,1/4,1/4,1] 
error rate 

20 

6.30e-01 


6.20e-01 


7.20e-01 


3.90e-04 


40 

3.06e-01 

1.04 

2.58e-01 

1.26 

2.83e-01 

1.34 

1.56e-04 

1.32 

80 

5.46e-02 

2.49 

5.93e-02 

2.12 

5.21e-02 

2.44 

4.27e-05 

1.87 

160 

7.52e-03 

2.86 

9.00e-03 

2.72 

8.37e-03 

2.63 

1.09e-05 

1.97 

320 

9.78e-04 

2.94 

1.18e-03 

2.92 

1.34e-03 

2.63 

2.74e-06 

1.99 

640 

1.23e-04 

2.99 

1.47e-04 

3.00 

4.03e-04 

1.73 

6.86e-07 

2.00 

1280 

1.54e-05 

3.00 

1.83e-05 

3.00 

9.11e-05 

2.14 

1.71e-07 

2.00 

2560 

1.92e-06 

3.00 

2.27e-06 

3.00 

2.34e-05 

1.95 

4.29e-08 

2.00 


1-norm error on the linear transport at t = 1 


N 

uniform 

error rate 

quasi¬ 

error 

unif 

rate 

random 

error rate 

[1,1/2,1/4,1/4,1] 
error rate 

20 

8.20e-02 


9.69e-02 


8.31e-02 


4.10e-02 


40 

2.75e-02 

1.57 

4.20e-02 

1.20 

2.76e-02 

1.58 

8.33e-03 

2.30 

80 

4.95e-03 

2.48 

9.47e-03 

2.15 

5.01e-03 

2.46 

1.25e-03 

2.73 

160 

7.35e-04 

2.75 

1.56e-03 

2.60 

7.40e-04 

2.75 

1.61e-04 

2.96 

320 

9.36e-05 

2.97 

2.10e-04 

2.88 

9.40e-05 

2.97 

1.94e-05 

3.06 

640 

1.14e-05 

3.04 

2.57e-05 

3.03 

1.14e-05 

3.03 

2.38e-06 

3.03 

1280 

1.41e-06 

3.01 

3.16e-06 

3.02 

1.41e-06 

3.01 

2.97e-07 

3.00 

2560 

1.76e-07 

3.00 

3.94e-07 

3.00 

1.76e-07 

3.00 

3.71e-08 

3.00 


Table 5: Discrete maximum norm error on numerical derivative (top) and 
discrete 1-norm error in linear transport equation with WEN03. e = h 2 , 
uq(x) = sin (2wx — sin(27r:r)/27r). 




Figure 3: Linear transport of smooth data on random grids. Discrete 1-norm 
error (left) and discrete maximum norm (right). 
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Maximum norm error on the numerical derivative on uq(x) 


N 

uniform 

error rate 

quasi¬ 

error 

unif 

rate 

random 

error rate 

[1,1/2,1/4,1/4,1] 

error rate 

20 

4.83e-01 


8.77e-01 


6.55e-01 


5.13e-04 


40 

3.14e-01 

0.62 

2.62e-01 

1.74 

3.06e-01 

1.09 

1.64e-04 

1.64 

80 

5.05e-02 

2.64 

5.48e-02 

2.26 

7.18e-02 

2.09 

4.32e-05 

1.92 

160 

5.46e-03 

3.21 

6.59e-03 

3.06 

8.61e-03 

3.06 

1.09e-05 

1.98 

320 

5.98e-04 

3.19 

7.58e-04 

3.12 

1.51e-03 

2.51 

2.74e-06 

1.99 

640 

7.10e-05 

3.07 

8.83e-05 

3.10 

2.77e-04 

2.44 

6.86e-07 

2.00 

1280 

8.73e-06 

3.02 

1.06e-05 

3.05 

6.43e-05 

2.10 

1.72e-07 

2.00 

2560 

1.09e-06 

3.01 

1.31e-06 

3.02 

1.61e-05 

1.99 

4.29e-08 

2.00 


1-norm error on the linear transport at t = 1 


N 

uniform 

error rate 

quasi¬ 

error 

unif 

rate 

random 

error rate 

[1,1/2,1/4,1/4,1] 
error rate 

20 

8.22e-02 


9.91e-02 


8.27e-02 


3.90e-02 


40 

2.40e-02 

1.78 

4.02e-02 

1.30 

2.41e-02 

1.77 

6.55e-03 

2.57 

80 

3.57e-03 

2.75 

7.68e-03 

2.38 

3.60e-03 

2.74 

8.54e-04 

2.94 

160 

4.57e-04 

2.97 

1.05e-03 

2.87 

4.64e-04 

2.95 

9.91e-05 

3.11 

320 

5.36e-05 

3.09 

1.25e-04 

3.06 

5.45e-05 

3.09 

1.07e-05 

3.20 

640 

6.35e-06 

3.08 

1.45e-05 

3.10 

6.49e-06 

3.07 

1.25e-06 

3.11 

1280 

7.80e-07 

3.02 

1.76e-06 

3.04 

7.97e-07 

3.02 

1.52e-07 

3.03 

2560 

9.72e-08 

3.00 

2.18e-07 

3.01 

9.94e-08 

3.00 

1.89e-08 

3.01 


Table 6: Discrete maximum norm error on numerical derivative (top) and 
discrete 1-norm error in linear transport equation with CWEN03. e = h , 
uq(x) = sin ( 2 nx — sin(27ra)/27r). 


24 







Figure 4: Linear transport of discontinuous data on random grids. Zoom of 
the solution close to the discontinuity (top), the discrete 1-norm and the total 
variation, respectively (bottom). 
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WEN03 - cf 1=1.55 



CWEN03 - cfl=1.25 



Figure 5: Stability region of SSP-RK3 and spectrum of first derivative operator 
computed with upwind flux and WEN03 reconstructions (left) and CWEN03 
reconstruction (right). 


puted the 1-norm of the error and the total variation: the results are reported 
in Figure [4j The errors in 1-norm are much closer to each other than in the 
smooth case, with e = h still being slightly better than the other choices. The 
test on the total variation shows that the increased resolution in the smooth 
case is obtained when the reconstructions stay closer to the central one and in 
fact the choices giving the lower errors in the previous test (i.e. CWEN03 and 
e = h) produce more total variations than the other ones. In any case the total 
variation is diminishing under grid refinement. 

4.3 Stability 

In Figure [5] we compare the stability region of the third order SSPRK used for 
time advancement and the spectrum of the operator that computes numerically 
the first order derivative in the approximation of the linear transport equation 
u t + u x = 0. The spectrum represented in the figure is the spectrum of the 
linearization of the nonlinear operator in the Fourier basis, which represents 
an analogue of the classical Von Neumann analysis that is applied to linear 
schemes. In particular, on a grid of 65 cells, we compute column-wise a 65 x 65 
matrix M as follows: for each column, we set up initial data coinciding with 
one of the real-valued Fourier basis with frequency between 0 and 32, compute 
the boundary extrapolated data, the numerical fluxes, the vector of differences 
u ~jj r \j 2 — an d decompose this latter quantity along the real-valued Fourier 

basis. The eigenvalues of M are shown in the picture, for WEN03 (left) and 
CWEN03 (right) and different choices of e. Linear stability is achieved if all 
the eigenvalues stay within the stability region of the Runge-Kutta scheme, 
i.e. in the white region of the plot. By symmetry, only the positive imaginary 
half-plane is shown. The CFL number was chosen so that the eigenvalues are 
very close to the boundary of the absolute stability region, showing that, for 
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average N cells average N cells 


Figure 6: Efficiency diagrams for u t + u x = 0 with periodic boundary conditions 
and smooth data. The initial datum is slowly varying on the left and quite 
oscillating on the right. 


both WEN03 and CWEN03, the choice e = h has a slight stability advantage 
over the other three. We also observe that the spectrum obtained with the 
CWEN03 reconstruction is less elongated in the imaginary axis direction than 
the one obtained with WEN03, as indicated by the different location at which 
the spectrum touches the boundary of the stability region of the Runge-Kutta 
scheme. 

4.4 Nonlinear conservation and balance laws 

Interplay with h-adaptive schemes The use of non uniform meshes is very 
important for non-linear conservation laws, since in the area around shocks, the 
truncation error can not be better than 0(h ) irrespectively of the order of the 
scheme employed and thus grid refinement is the only available tool to reduce 
the computational error in that area. The results of this paper are relevant for 
moving mesh methods (quasi-uniform grids) and for locally refined a/lyd-grids. 
Here we present tests in the latter setting, using the third order generalization 
of the h-adaptive scheme presented in [21], which was analized in [23 . The 
scheme employs a single non-uniform mesh, built from an initial uniform one 
by locally and recursively splitting in two the troubled cells according to some 
error indicator. As in m, the numerical entropy production is employed as an 
error indicator, but we believe that this result are rather independent on the 
details of the adaptive strategy, since in any case the numerical scheme would 
have to deal with nearby cells whose size ratio is a power of two and are thus 
not quasi-uniform. 

In Figure [6] we present the results obtained integrating the linear transport 
equation with periodic boundary conditions on [0,1] and initial datum u\[x) = 
sin(7ra; — sin(7rx)/7r) (left) and U 2 {x) = sin(7ra;) + 0.25sin(157ra;)e _20a: (right). 
The grid is adaptively refined during the numerical integration, using cells of 
three (four) sizes respectively, thus with maximum grid ratio of 1/8 or 1/16. 
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Figure 7: Efficiency diagrams for the Burgers’ equation. The initial datum on 
the left is M 3 and gives rise to a standing shock (final time 0.35), while the one 
on the right is W 4 , which is much more oscillating and gives rise to a complex 
solution structure with moving shocks (final time 0.45). 


The discrete 1-norm error is plotted against the average number of cells that 
were used during the computation, which is a proxy of the computational effort 
of the scheme. Since the flow is smooth, in this test there are no theoretical 
reasons why non-uniform grids would outperform uniform ones. However, the 
figure shows that choosing an h-dependent e yields consistently lower errors and 
that the choice e = h clearly outperforms the other two. 

In Figure[7]we present the results obtained integrating the Burgers’ equation 
ut + u 2 ) x = 0 on the domain [0,1] with the initial datum 1 * 3 ( 2 ;) = — sin( 7 T 2 ;) 
(left) and u^(x) = — sin( 7 T 2 :) + 0.2sin(5.07r:r) (right). The adaptive scheme 
clearly converges at a faster rate than the uniform grid one which is locked 
into first order behaviour by the presence of the shocks. In the tests shown on 
the left, 113 gives rise to a standing shock in the middle of the domain and we 
employed a coarse grid of 16-2 fc cells with 3+2A; cell sizes (k=0... 3) as this is the 
choice that could recover third order convergence with respect to the average 
number of cells (see [25] 1. In the test on the right, U 4 gives rise to a richer 
solution structure with two moving shocks converging towards the middle and 
four little waves where third order accuracy plays an important role. Here we 
employed 16 • 2 k cells with 4 + 2 k cell sizes (k=0,..., 4). In both cases, choosing 
an h-dependent e yields consistently lower errors and the choice e = h clearly 
outperforms the other two, especially in the more complex situation depicted 
on the right. 

Finally, we consider the classic problem from [25] of the interaction of a shock 
with a standing acoustic wave. The conservation law is the one-dimensional 
Euler equations with 7 = 1.4 and the initial data is 

f (3.857143,2.629369,10.333333), x G [0,0.25], 

(P,+.P) - j ( L0 + 0 . 2 sin( 167 rx), 0 . 0 , 1 . 0 ), x G (0.25,1.0], 

The evolution was computed up to t = 0.2 and the results are presented in 
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10 2 10 3 10 4 
average N cells 


Reference solution 



x 



Figure 8: Shock-acoustic interaction test. Top-left: efficiency diagram. Top- 
right: plot of the density (reference solution). Bottom: comparison of the 
numerical solutions obtained with different choices of e in the areas marked 
with dashed rectangles in the top-right panel. 
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N = 100 

N = 200 

N = 400 

N = 800 

WEN03+P2 

CWEN03 

1.63 • HE 2 

6.29 • 10~ 4 

1.09- 10" 2 
1.27- 10" 4 

1.51 • 10~ 2 
3.00 • HT 5 

1.06 • 10” 2 
6.88 • 10~ 6 


Table 7: Test on the violation of the passivity constraint. Conservation errors of 
numerical schemes modified by truncating to 0 all negative water height values. 


Figure [8j As the right moving shock impinges in the stationary wave, a very 
complicated smooth structure emerges and then gives rise to small shocks and 
rarefactions (top right). The evolution was computed with uniform grids and 
also with the adaptive algorithm, using coarse grids of 32,64,128 cells, with 
respectively 6, 8,10 refinement levels. The plot of the error versus the average 
number of cells employed in each run shows the effectiveness of the h-adaptive 
procedure and the considerably lower errors obtained when e = h is employed 
in the reconstruction procedure (top-left). Three significant portions of the 
numerical solutions obtained with the h-adaptive algorithm (128 coarse cells 
and 10 levels) are shown in the lower part of the figure. It is evident that 
the choice e = h gives the best results in the approximation of both the high 
frequency waves and of the little shocks on the left, in particular avoiding the 
formation of spurious waves, which occour when using e = 10~ 6 or e = h 2 . 

Third order well-balancing We now consider the finite volume discretiza¬ 
tion of the shallow water equations. It is well known that a second order ac¬ 
curate scheme which is well-balanced for the lake at rest steady state can be 
constructed with the classical hydrostatic reconstruction technique of [2|. In 
this scheme, the cell average of the source term is computed as 

= \ (kj-l/2 + h j + 1 / 2 ) (z j+ 1/2 - Zj_ 1 / 2 ) 

where hj± 1/2 and Zj± 1/2 are the point values for the water height and the bottom 
elevation at the boundaries of the j-th cell that were employed in the compu¬ 
tation of the numerical fluxes. A third order accurate scheme may be obtained 
replacing the original linear, slope-limited, reconstruction with an higher order 
non-oscillatory one and replacing the original second-order accurate quadrature 
of the source term with a third-order accurate one. A simple way to obtain such 
a formula is to apply the Richardson extrapolation technique to Si, as proposed 
hi |20| . Namely, denoting with S 2 the same quadrature rule applied with two 
sub-intervals, we have that the combination ( 4 S 2 — Si)/3 is a fourth-order ac¬ 
curate quadrature for the source term. However, since 

£2=5 (hj~ 1/2 + hj) (Zj - zj_ 1/2) + \ (hj + h j+ 1/2) (~j+i/2 - Zj) 

the point-value reconstructions at the cell center must be available as well. 

As mentioned at the beginning of (j3j WEN03 cannot be applied to com¬ 
pute third order accurate reconstruction at cell centres hj and Zj. Since these 
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values are needed in the well-balanced quadratures, one could use hj = hj and 
Zj = Zj instead of the third order accurate value, but this would reduce the 
convergence order to 2 as it can be easily checked on the smooth flow sug¬ 
gested in Alternatively, one may reconstruct the point values hj and Zj 
by evaluating the central optimal polynomial at the cell center. Let’s call this 
reconstruction WEN03+P2 and compare it with CWEN03. Neither CWEN03 
nor WEN03+P2 give rise to positive schemes, but the latter causes more os¬ 
cillations, expecially close to wet-dry fronts. In order to show the differences, 
we consider and compare the positive schemes obtained by artificially setting 
/i” +1 = 0 whenever /i” +1 < 0 after each timestep. These schemes are of course 
not conservative, but their conservation error accumulates over time and gives 
a measure of the violation of the positivity constraint. Table [7] shows the con¬ 
servation errors registered for the simulation of a pond with bottom z(x) = 2a; 2 , 
initial water level H(x) = z(x) + h{x) = max(l + 0Ax,z(x)) and q(x) = 0. 
At time t = 4.0 the water surface has inverted its movement four times, each 
time giving rise to spurious waves at the dry points. The scheme based on 
WEN03+P2 does not show any convergence of the conservation error, indicat¬ 
ing that the positivity constraint violations are much severe than those caused 
by the CWEN03 scheme, which gives rise to a conservation error that decays 
as N~ 2 with the increase of the number of points. The results in the Table [7] 
refer to the case of quasi-uniform meshes, but we point out that similar ones 
were observed on uniform ones. 


4.5 Two-dimensional tests 


A generalisation of the CWENO reconstruction to locally-refined quad-tree grids 
was described in [23j . Here we present a comparison of the different choices for 
e in this setting. 

The reconstruction presented in (231 is a truly two-dimensional reconstruc¬ 
tion that generalises the one of |T5] to non-uniform meshes. In brief, for the cell 
ttj one considers as “optimal” second order polynomial P J OPT (x, y) that interpo¬ 
lates the cell average in Qj exactly and all the cell averages in the neighbours in a 
least-squares sense (this approach is of course needed since in a two-dimensional 
h-adapted grid the number of neighbours is variable). Additionally the recon¬ 
struction considers four “directional planes”, i.e. first order polynomials that fit 
exactly the cell average in and, in a least squares sense, the cell averages of 
the neighbours in the north-east, north-west, south-west and south-east sector. 
For the precise definition of the neighbourhoods and stencils, see [25J. Then a 
second order polynomial P^(x,y) is defined by requiring that 


P° PT (x,y) = Ijf (x,y) + hP^(x,y) + lPr W (x, V )+ lPf W (x,y)+ iP^fay) 


non-linear weights are computed from the optimal weights Cc = \ and Cne = 

Cnw = Csw = Cse = | with the help of the Jiang-Shu indicators and a 
reconstruction polynomial is defined by 

P REC (a;, y) = wc Pf (x, y)+uJNEP™ E (x, y)+u> NwPf W (x, y)+^swPf W (x, y)+wsEPf E (x, y). 
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Inin 

SWSii 


1 

N 

e(h) = 
error 

Ttr 5 

rate 

e(h) = 
error 

h 2 

rate 

e(h) = 
error 

h 

rate 

Qo 

1.03e4 

1.82e-02 


1.30e-02 


2.94e-03 


Q i 

4.13e4 

4.59e-03 

1.98 

1.93e-03 

2.75 

2.03e-04 

3.86 

02 

1.65e5 

1.10e-03 

2.06 

2.16e-04 

3.16 

2.31e-05 

3.13 

03 

6.62e5 

1.20e-04 

3.19 

2.04e-05 

3.40 

2.91e-06 

2.99 

04 

2.65e6 

4.97e-06 

4.59 

2.08e-06 

3.30 

3.65e-07 

3.00 


Table 8: Maximum norm errors for the two-dimensional CWENO reconstruction 
on h-adapted grids. The error decay rates are computed with respect to y/N, 
where N is the number of cells in the grid. (Due to the choice of grids Qk, N 
scales as the inverse of the size of the smallest cell in the grid). The grid Qq and 
the function u is depicted on top. 


This latter polynomial is uniformly third order accurate in the cell (for smooth 
data) and it can be evaluated where needed. For non-smooth data, the pro¬ 
cedure of computing the nonlinear weights from the linear one by using the 
regularity indicators, ensures that data from non-smooth sub-stencils are not 
relevant for the coefficients of P^ 0 and the reconstruction is non-oscillatory. 
For more tests, see [23). 

In this paper we want to test the influence of the choice of e in the accuracy 
of the reconstruction for smooth data. We remark that the role of e is quite 
important in h-adapted meshes, since the grid includes cells of very different 
size and it is thus quite difficult to choose a fixed value of e that works well on 
every cell. 

To this end we consider the function u(x, y) = sin(27rcc) cos(27 xy) on the 
unit square with periodic boundary conditions on a uniform coarse grid. The 
grid is locally refined with a simple gradient indicator, obtaining the locally 
adapted grid Q 0 depicted in Table [8] Finer grids Qk are then obtained from 
this one by splitting each cell of Qo into 4 fc equal parts. For each grid, the 
cell averages of u are computed with a fifth-order gaussian quadrature rule and 
then the reconstruction of boundary extrapolated data is performed using the 
two-dimensional CWENO reconstruction. The reconstruction polynomials are 
computed in each cell and evaluated at the quadrature nodes of the third-order 
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gaussian quadrature, which would be the quadrature of choice for computing 
fluxes of a finite-volume scheme on such grids. These values are then compared 
with the exact values of the function u and the maximum-norm error is reported 
in Table [8] 

We observe that the choice e(h) = h gives the best results both in terms of 
error decay rates and in terms of absolute values of the error. In fact, the last 
column of Table [8] consistently records lower errors than the other two and has 
a more regular behaviour of the decay rates. 

5 Conclusions 

HJ and Q3] discussed the advantage of choosing the parameter e in WENO 
and, respectively, CWENO boundary value reconstruction procedures. In this 
work we extended their results to non-uniform meshes and concentrate on finite 
volume schemes, whereas [TJ worked in the finite differences approach. 

Our work shows that, also for a non-uniform grid, choosing e as a function 
of the local mesh size (i.e. e = e(hj) for the reconstruction of u+_ x , 2 and 
U J+ 1 / 2 ) all° ws one t° recover the optimal error of convergence even close to 
local extrema of the function being reconstructed and in general provides a 
much more regular pattern of error decay, which is a quite important feature 
for the good performance of error indicators. 

We compared the choices e = h and e = h 2 , showing that both can yield the 
aforementioned improvements and that the former is slightly better on smooth 
tests, while the latter gives a slightly lower total variation increase in discontin¬ 
uous cases. 

Furthermore, we compared the performance of the different choices for e in 
an h-adaptive context, with third order scheme that employs the CWEN03 
reconstruction and, at each timestep, locally refines and coarsens the grid ac¬ 
cording to the numerical entropy production error indicator. We observed that 
reconstructing in each cell flj with ej = hj (where hj here is the local cell size) 
yields the best results also in the shocked cases. The paper is completed by a 
test on the performance of the CWEN03 reconstruction for the shallow water 
equation with dry states. 

Finally, we point out that €j = hj has been already successfully employed 
in the h-adaptive scheme for conservation laws in two space dimensions based 
on a two-dimensional generalisation of the CWEN03 reconstruction on locally 
adapted meshes presented in [23] and that more testing on the comparison of 
different choices for e is presented here. Since that reconstruction is based on 
local approximation and not on interpolation, the analysis of that more general 
situation goes beyond the scope of this paper, but constitutes an interesting line 
of research. 
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